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Abstract 

We present an energy expression for restricted open-shell Kohn-Sham theory for N unpaired 
electrons and single-electron operators for all multiplets formed from up to five unpaired electrons. 
It is shown that it is possible to derive an explicit energy expression for all low-spin multiplets of 
systems that exhibit neither radial nor cylindrical symmetry. 



I. INTRODUCTION 



Restricted open-shell theory [H [2] and in particular restricted open-shell Kohn-Sham 
theory (ROKS) [21 HI |S] has gained renewed interest in recent years due to its application in 
the simulation of photoreactions [6l [71 El [9l HOI E] • While the concept of restricted open- 
shell theory is seemingly simple, only part of the states of interest are accessible in present 
implementations. For the high-spin case (multiplicity = number of unpaired electrons + 
1) the situation is simple: only one Slater determinant is needed for the description of the 
wavefunction. The most trivial case is the one with but one unpaired electron which gives 
rise to a doublet that is described by a single determinant in a straightforward way. For 
most single-determinant cases, results very similar to the unrestricted theory are obtained 
and the use of the restricted theory is not necessary. For low-spin cases the situation 
is much more involved and from a purely mathematical view one might doubt if these 
cases can be analysed at all. One of the general criticisms concerning open-shell theory 
has already been addressed by Roothaan in his original restricted open-shell Hartree-Fock 
(ROHF) paper \T\: He showed that although orbital-dependent operators are obtained, 
the equations can be augmented to a Hermitean formulation by the use of projection 
operators. This guarantees diagonalization to real eigenvalues. For low-spin cases there 
is the additional problem of orbital rotations leading to unphysical localization which has 
also been noted as sudden polarization in related approaches |T2j, and there is a lot of 
confusion in the literature due to the difficulty in correctly interpreting and avoiding this 
phenomenon of unphysical orbital localization. Undesired rotations can in be avoided by 
exploiting spatial symmetry or by state averaging; general self-consistent solutions for the 
restricted open-shell problem are difficult to achieve [T3]. In recent years we have shown 
for the case of restricted open-shell Kohn-Sham theory, that by proper modification of the 
off-diagonal elements of the Kohn-Sham matrix these localizations can even be avoided in 
self-consistent calculations for first-principles molecular dynamics simulations [H [S]. The 
parameters used in these algorithms are case dependent (localized / delocalized situations); 
up to now no unique approach exists that would efficiently converge to the correct solution 
for all cases. Finally, if one wants to use Kohn-Sham instead of Hartree-Fock expressions, 
there is the question how the energy expression for a multi-determinant approach should be 
determined. In ROKS [3] we use the energy expression as proposed in the sum method [H] . 
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This expression reduces to the ROHF terms if it is used with the exact energy expectation 
values of the open-shell Slater determinants, instead of inserting the Kohn-Sham expressions. 

Practical ROKS dynamics calculations are presently restricted to the case of two unpaired 
electrons which is by far the most important case for the description of photoreactions. This 
restriction has several reasons: 1. While the concept is clear, there is, to our knowledge, no 
explicit formulation of the energy expression for arbitrary low spin states in the Kohn-Sham 
literature. 2. Based on an energy expression, single-electron ROKS equations have to be 
derived for the particular spin densities involved. 3. Algorithms for the self-consistent 
solution of the equations for N electrons must be developed and tested. In the present 
paper, we address points 1 and 2. 

We simplify the derivation, which corresponds to deriving an energy expression for a 
spin-adapted configuration, by neglecting the possibility of symmetry-determined degen- 
eracy as it may be present in atoms = 0) or in diatomic molecules {[H^Lz] = 0). 
The treatment of such highly symmetric systems hardly plays a role in molecular dynamics 
simulations. Instead we use the occupation pattern as a symmetry like suggested in a 
similar way by the work of Ziegler, Rank and Baerends p!^, Daul [15], and Noodleman 
[T6| IT7j . On this basis it is possible to compute, for example not only the lowest singlet 
state, but the singlet states with zero, two, and so on, unpaired electrons. (Note that the 
singlet state with zero unpaired electrons is not necessarily the lowest one.) In calculations 
based on the sum method, in addition to the occupation pattern high spatial symmetry was 
exploited for the computation of multiplets in the context of ligand field theory (see [18] 
and references cited therein). We do not make use of radial or cylindrical symmetry (or of 
approximate radial symmetry, if a central metal atom is in an environment). This means 
that for example, for two unpaired p electrons we cannot compute ^D, ^D, ^P, ^P, ^S, and 
-^S, but just Ti and Si. Higher states are only accessible as local minima. Ignoring the fact 
that an atom might be in an environment with a particular symmetry seems appropriate 
for reactive molecular dynamics simulations, where the spatial symmetry and the degree of 
degeneracy should be free to change at any time. Note that the restriction of not treating 
symmetry-determined degeneracy explicitly, does not mean that degeneracy cannot be 
described at all. It just means that degeneracy is determined at the single-configuration 
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level - which may turn out to be wrong in comparison with a multi-configuration approach. 

In the next section, the energy expression for N unpaired electrons is derived. In the 
following two sections, first the general form of the ROKS operators for N electrons is given 
and then the explicit symmetry determined parameters of these operators are specified for 
up to five unpaired electrons. 

II. ENERGY EXPRESSION FOR N ELECTRONS 

SD C 




FIG. 1: Sketch of the relations between single-determinant levels and single-configuration levels 
for an N electron case. The energy levels for a certain multiplet i (right side) are 
degenerate, the single-determinant energies are averaged for a certain j (corresponding 
to a certain number of a/P electrons, left side). Hence the indices I and k numbering 
the different permutations / micro states can be omitted in the derivation of the general 
energy expression. 
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Since [H, S"^] = 0, an open-shell electronic configuration "^f must in general be composed 



of several Slater determinants "^j^ : 



In this way it is possible to describe a spin-pure energy state. For the derivation of the 
energy expression, the coefficients dij do not have to be explicitly determined. 
The energy expectation values of a configuration \E'p with N unpaired electrons can be 
obtained as a sum of the energy expectation values of the single determinants: 

Ef-" = j:c,Ef-'' (2) 

j 

In the derivation of the coefficients Cij we use the short notation: 

E^ = j:c,Ef- (3) 

j 

with 



2 

where n{a) and n{(3) are the total numbers of a and (3 electrons respectively. 
The Slater determinants which can be formed for a particular number of unpaired electrons 
N, can be ordered according to their magnetic spin quantum number Ms: 

Ms = linia)-nm (5) 
With this ordering the Hamilton matrix for the wavefunction [T] factorizes (Fig. [2]). 
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FIG. 2: Hamilton matrix, ordered according to the Ms values of the Slater determinants. Only 
the elements in the colored blocks differ from zero. 

For Smax and -S^ax (high-spin: all unpaired electrons have either a oi (3 spin) the 1x1 
block consists of the expression for only one Slater determinant and = E^^q = EfP^^g^^^ 
for the high-spin case. From this, an energy expression can be obtained for the block of the 
matrix with Ms = Smax — 1: 
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SD \ 

M S=Smax — 1 / 
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S — Smax 1 



+ 



N 



S — Smax 



(6) 



Similarly, for the block of the matrix with Ms = Smax ~ 2: 



2 ; {^Ms=s„ 



pS—Smax —2 



(7) 



With equation [6] it follows: 



2 l\-^Ms = Smaa^-2/ 



E^ 



S — Smax 2 



+ 1-1 (^Ms=Sma.-l) 



(8) 



This result can be generalized for the matrix block with Ms = Smax — j'- 



SD \ 

Ms = Smax—j I 



N 

i-i 
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or: 



■s=s, 



(10) 



In a shorter notation, we write: 



E. 



c 



(11) 



From this formula it follows: 



E. 



C ^N+l^ SD^ 

^ N + l-2r ^ ' 



3 I TpSD \ 

Ar + i_2/Vi) 



(12) 



That is, the energy of a state with multiplicity M = 2j + 1 can be determined from the 
energies E^^ and E^P^, all other coefficients Cj are zero. 
The resulting energy levels are depicted true to scale in Fig. |3| 
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FIG. 3: Energy diagram for up to five unpaired electrons. In the determination of the energy 



levels, equation 12 was used. The energy levels are numbered by the index j. 
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For the computation of multiplet energies, Noodlemans elegant and simple formula is 
frequently used [T5| : 



= (1 + c)E^^ - cE^^ (13) 

with 

' = s Is 

In our notation this equation reads: 

E^ = {1 + c)Ef - cEl'' (15) 

In contrast to equation 12, the energies are determined from E^^ and Eq^ . The Noodle- 
man formula results in identical energy terms for up to three unpaired electrons. (Two 
unpaired electrons: triplet: Eq = Eq^, singlet: E^' = 2Ef^ — E^^, three unpaired elec- 
trons: quartet: Eq = E^^ , doublet: = \Ef^ — ^E^^). For more than three unpaired 
electrons, as is readily verified, the Noodleman formula does not fulfill the sum rule which 
says that the sum of the energies of the configurations must equal the sum of the energies 
of the determinants they are formed from. 

III. GENERAL ROKS EQUATIONS 

Starting from 

Ef = c,E, + c,^,E,_, (16) 
we derive single-electron ROKS equations. Using the coefficients Cj as determined by 



equation 12 , the energy can be written in terms of energy expressions of single Slater deter- 
minants, whereby summation over half the determinants is sufficient since there are always 
two determinants with equal energy for which the a and /5 spins are just interchanged. 
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Ef =c, 



2 



(17) 



+ c 
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Using the Kohn-Sham energy expression 
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(18) 
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the general ROKS operators for shell i are obtained by functional variation: 



2 



+ 2^ 



1 (n-,^;-Jp-,pfj+nf^,,;^,[p-,pfj] 



Jl 

2 



(19) 



3-1 

2 



More details of the derivation are given for the case of two unpaired electrons in 



For practical applications of equation 19, the Uj^; have to be determined, whereby the index 
i denotes the - closed and open - shells, the index j which numbers the energy levels is 
defined in equation |4| the index / numbers the microstates belonging to a certain j. The 
coefficients njjj are given in the following section for up to five unpaired electrons. 



IV. SPECIAL ROKS OPERATORS 



A. Two unpaired electrons 



The case of two unpaired electrons is described in detail in [3]. The energy expressions 
are as follows: 
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Eq 



(20) 

= 2Ef^ - E^"" (21) 

The spin densities are given in Fig. |4} The prefactors in equation 19 are hsted in table [T] 
for the non-trivial singlet case. 

A B 



A 



la lb la lb 

FIG. 4: Spin densities for the determinants that can be formed for two unpaired electrons, 
grouped according to their value of j. Here and in the following figures, the indices 
1 are described with la, lb, 2a, 2b etc., indicating that there are always two determi- 
nants with the same energy. 
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TABLE I: Coefficients in equation 
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for two unpaired electrons. 



This describes Kohn-Sham operators which explicitly read: 
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Triplet case: 

Pol 

Singlet case: 
F, 
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B. Three unpaired electrons 
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FIG. 5: Spin densities for three unpaired electrons. 
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TABLE II: Coefficients in equation 



19 



for three unpaired electrons. 



Tlie Kolin-Sliam operators for tliree unpaired electrons (quartet and doublet) are: 
Quartet case: 



= - -V^ + 
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(26) 
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Doublet case: 
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(27) 



^^xApIii Pai\ 

To simplify the representation for the cases with more electrons we introduce a short 



notation for these equations (Tables III and IV). 
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TABLE III: ROKS operators for three unpaired electrons, quartet case. 
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TABLE IV: ROKS operators for three unpaired electrons, doublet case. 



Four unpaired electrons 
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FIG. 6: Spin densities for four unpaired electrons. 
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Table [V] lists the coefficients for four unpaired electrons, tables VI, VII and VIII the Fock 
operators in short notation. 
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TABLE V: Coefficients in equation 
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for four unpaired electrons. 
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TABLE VI: ROKS operators for four unpaired electrons, quintet case. 
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TABLE VII: ROKS operators for four unpaired electrons, triplet case. 



15 











c 
















B 














1 


2 


3 






1 


2 


3 


4 






in 


+ f[|( 


a 




a 


/3 


a 


/3 


)]- 




a 


/3 


a 


/3 


a 


13 


Q 


/3 


)] 


c 


2 




1 


1 


1 


1 


1 


1 






1 


1 


1 


1 


1 


1 


1 


1 




ol 


1 




1 





1 





1 









1 





1 





1 








1 




o2 


1 




1 








1 





1 






1 





1 








1 


1 







o3 


1 







1 





1 


1 









1 








1 


1 





1 







o4 


1 







1 


1 








1 









1 


1 





1 





1 








TABLE VIII: ROKS operators for four unpaired electrons, singlet case. 



D. Five unpaired electrons 
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FIG. 7: Spin densities for five unpaired electrons. 



Table IX lists the coefficients for five unpaired electrons, tables |X, XI and XII the Fock 
operators in short notation. 
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TABLE IX: Coefficients in equation 
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for five unpaired electrons. 
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TABLE X: ROKS operators for five unpaired electrons, sextet case. 
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TABLE XI: ROKS operators for five unpaired electrons, quartet case. 
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TABLE XII: ROKS operators for five unpaired electrons, doublet case. 



V. CONCLUSIONS 

We have derived a general explicit energy expression for restricted open-shell Kohn-Sham 
theory which fulfills the sum rule for the energies. If degeneracy is not prescribed and the 
occupation pattern is interpreted as a symmetry, the energy expression is also valid for the 
sum method by Ziegler, Rank, and Baerends, and, if the exact energy expectation values are 
used, for ROHF itself. In the latter case the energy expression reduces to the exact energy 
of a single configuration [T9] . 

By inserting the Kohn-Sham expressions, we have derived ROKS operators and have given 
explicit expressions for these operators for up to five electrons. 

The energy expression and the operators constitute just part of the restricted open-shell 
problem. For more than one open shell, only the high-spin solution is easily obtained. For 
the low-spin solutions of the ROKS equations, specific SCF algorithms turned out to be 
useful in the case of two open shells. In future work we want to extend these algorithms to 
more open shells. 
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